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ABSTRACT 

We present properties of moderately massive clusters of galaxies detected by the newly developed Hyper 
Suprime-Cam on the Subaru telescope using weak gravitational lensing. Eight peaks exceeding a S/N ratio of 
4.5 are identified on the convergence S/N map of a 2.3 deg^ held observed during the early commissioning 
phase of the camera. Multi-color photometric data is used to generate optically selected clusters using the 
CAMIRA algorithm. The optical cluster positions were correlated with the peak positions from the convergence 
map. All eight signihcant peaks have optical counterparts. The velocity dispersion of clusters are evaluated by 
adopting the Singular Isothemal Sphere (SIS) ht to the tangential shear prohles, yielding virial mass estimates, 
MsoOc^ of the clusters which range from 2.7 x 10*^ to 4.4x 10'^ Mq. The number of peaks is considerably larger 
than the average number expected from ACDM cosmology but this is not extremely unlikely if one t^es the 
large sample variance in the small held into account. We could, however, safely argue that the peak count 
strongly favours the recent Planck result suggesting high erg value of 0.83. The ratio of stellar mass to the dark 
matter halo mass shows a clear decline as the halo mass increases. If the gas mass fraction, fg, in halos is 
universal, as has been suggested in the literature, the observed baryon mass in stars and gas shows a possible 
deheit compared with the total baryon density estimated from the baryon oscillation peaks in anisotropy of the 
cosmic microwave background. 

Subject headings: cosmology: observations—dark matter—large scale—structure of universe 


1. INTRODUCTION 

Clusters of galaxies are the largest gravitationally bound 
systems in the universe and have been useful cosmological 
probes to learn about the geometry and structure of the uni¬ 
verse. X-ray observations have, to date, provided the most 
efficient way to collect samples of clusters. Cluster cata¬ 
logs c ompiled from the ROSAT All-Sky Survey (IVoges et al.l 
1 19991) have often been used in combination with follow-up 
observations by XMM/Chandra to give rigorous cosmological 
constraints. Based on a sample of 238 clu sters with redshifts 
ranging from 0 to 0.5. lMantz et akl (120101) present constraints 
on the density parameters (Um and ug) from cluster abun- 
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dance and on the dark energy parameters (fix and wq) from 
the redshift evolution of that abundance. The mass fraction 
of host gas, /gas, measured at large radii in the most massive 
halos is also used to as a cosmological probe, as it provides a 
direct measure of the luminosity distance under the assump¬ 
tion of the universality of the gas fraction in massive halos; 
it is expected to scale as /gas(z) oc and is sensitive 

enough to give a level of constraint o n the parameters com- 
parable to those from other methods (lAllen et al.ll2008h . So 
far, no signihcant departure from ACDM model is reported. 
The eROSITA satellite is scheduled to launch in 2016 and is 
expected to collect yet more clusters in a more distant redshift 
range, and is expec ted to obtain more stringent constraints 
(iPredehl et al.ll2^ . 

In the mean time, the clusters themselves draw as¬ 
tronomers’ interests as a site of galaxy formation and galaxy- 
galaxy interactions. Galaxy formation is a process which is 
accompanied by the conversion of baryons into stars in a dark 
matter halo; the star formation efficiency is critically depen¬ 
dent on the ability of the gas to cool. In a less massive dark 
matter halo, the gas is easily reheated and swept out from the 
center by feedback processes such as supernova explosions. 
In a massive halo, the cooling time becomes longer as the halo 
mass becomes larger and the virial gas temperature higher, 
and galaxy formation is expected to be suppressed. It is also 
argued that AGN could provide an efficient star-formation 
suppression feedback mechanism in the higher mass range. 
These considerations indicate that there is an optimal halo 
mass where star formation is most effective around 10*^ Mq 
and where the stellar mass to halo mass ratio should peak. Sta¬ 
tistical estimate based on the abundance matching technique 
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gener ally agree with these naive expectations. (iBehroozi et alJ 

[ 2 ^ . 

In terms of individual halos, we now have accumulating 
observational evidences that the star formation efficiency ac¬ 
tually drops as the host dark matter halo mass increases 
in the group to cluster scale, > 10*^ Mr?,. (iLin et akl 
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iGonzalez et al.ll2013ll . These observa- 


2012 


tional efforts are also motivated by the desire to measure the 
total baryon mass (gas and stars) to dark matter ratio and to 
compare with the universal baryon to dark matter density ratio 
measured from the anisotropy of the cosmic microwave back¬ 
ground. This is, o f course, related to the interesting ’’mi ssing 
baryon problem” (iFukugita et al.lll998t lFukugitall2003l) . and 
quantitative measurements are important. 

However, the discrepancy among the observations are quite 
large so far. On the one hand, using 12 galaxy clusters at red- 
shifts around 0.1 with M = 1 ~ 5 x 10'^ IGonzalez et ^ 
(120131) argue that the difference between the universal value 
and cluster baryon fractions is less than the systematic uncer- 
tainti es associated with the m ass determinations. On the other 
hand, iLeauthaud et al.l (1201 2|) insists clear shortfall based on 
the observation of X-ray groups with 1 ~ 7 X 10'3 Mq, in 
COSMOS field. 


This confusion may arise par tly from the var i ety of halo 
mass measurement techniques. IGonzalez et al.l (120131) esti¬ 
mate the halo mass from the X-ray temperature using the 
usual yrrial scaling relation in which they assume the X- 
ray gas is in hydr odynamic equilib rium. The dynamical 
method adopted in lAndreonI (1201 Ol) assumes that the sys¬ 
tem is yrrialized, which many no t be tru e in the outer region 
of the clusters. ILeauthaud et al.l (120121) estimates the mass 
primarily from X-ray luminosity although this is calibrated 
by weak lensing. The chosen methods work best in differ¬ 
ent mass ranges, and their systematics almost certainly ex¬ 
ist and are almo s t certa inly different from technique to tech¬ 
nique. IWu et al.l (1201 5l) reported that the gas fraction is anti- 
cotTelated with stellar mass fraction in their simulation. This 
could shift the ratio of the stellar mass oyer the halo mass 
which introduces another complexity in the studies based on 
the optically or X-ray selected clusters. 

In this paper, we employ new approach to sample clusters 
and measure the halo mass directly by weak lensing with the 
fewest assumptions about the dynamical state of the cluster. 
As a part of the commissioning run on the S ubaru telescope 
of the newly deyeloped Hyper Suprime-Cam (iMiyazaki et al.l 
1201 2|) camera, we obseryed a 2.3 deg^ field in i-band. We 
use these obseryations to locate massiye dark matter halos di¬ 
rectly yia weak lensing by means of the deriyed conyergence 
map. A deep multi-color catalog is used to generate an op¬ 
tically selected cluster catalog with estimat ed redshifts using 
the noyel CAMIRA algorithm (IOgurill20T4h . This optical cat¬ 
alog is then correlated with the list of peaks in the conyergence 
map. These peaks, if not coincident with optically selected 
clusters, can be spurious, generated by noise, can correspond 
to chance coincidences of less massiye halos along the line 
of sight, or can corresponds to real clusters with yery high 
mass-to-light ratio. 

We chose the Deep Lens Suryey (DLS) Field 
dWittman et al.l l2002h for this inyestigation, for which a 
deep multi-color galaxy catalog is publicly ayailable, to 
generate the optical cluster catalog. Based on the shear 
selected cluster catalogs combined with the luminosity 


TABLE 1 

Exposure information 


Field 

Filter 

J2000 

Tcxp [sec] 

Med. Seeing [”] 

South 

HSC-i 

(139.50, 30.00) 

2700 

0.58 

North 

HSC-i 

(139.50, 30.75) 

3000 

0.57 


based catalog, we determine the cluster number count and 
the redshift distribution oyer a wide contiguous area oyer 
which obseryational systematic errors can be expected to be 
minimized, since the data are taken almost simultaneously 
oyer the whole field with an instrument with a wide field 
and relatiyely uniform Point Spread Function (PSF). We 
then determine the stellar mass fraction of indiyidual shear 
selected clusters. The indiyidual halo mass, estimated 
di rectly through lensing, rang es nicely between the samples 
of iLeauthaud et al.] (120121) and IGonzalez et all (l2013h . 

Not all the obserying facilities allow this approach. In or¬ 
der to locate the indiyidual dark matter halos, a high resolu¬ 
tion surface mass density map is necessary, which requires a 
sufficient number density of faint background galaxies whose 
shapes are measured with the necessary accuracy. With a 
much wider field of yie w (roughly ten times th e area) than the 
original Suprime-Cam (IMiyazaki et al.ll^02bl) . the 1.5 degree 
Hyper Suprime-Cam (HSC) field on the 8.2 m Subaru tele¬ 
scope has a crucial adyantage for a weak lensing suryey. We 
haye an approyed plan to ca rry out a three layer leg acy suryey 
using Hyper Suprime-Cam (IMiyazaki et al.ll^l3l) Based on 
the pilot obseryations presented here, we attempt to demon¬ 
strate the prospects of the power of the legacy suryey. 

2. DATA ANALYSIS 
2.1. Data set 

DLS Field 2 was obseryed on February 4th, 2013 Hawaii 
time. As is shown in Tabl^Il two pointing centers are chosen 
whose angular distance is 0.75 degree apart. Because the field 
of yiew of HSC is 1.5 degree, we haye substantial oyerlap be¬ 
tween two pointings. The exposure time is either 300 sec (9 
exposures on south, 8 exposures on north) or 150 sec (4 ex¬ 
posures on north), and the total exposure time of south and 
north pointing is 2700 sec and 3000 sec, respectiyely. A cir¬ 
cular dithering pattern of radius 2 arcmin is adopted around 
each pointing center to fill the gap of CCDs. The position 
angle is increased by 72^^90 degree between exposures. 



Fig. 1.— FWHM map (a) and ellipticity whisker map (b) of one of the 
exposure on the south (200 second). On the right, dash line shows the tilt axis 
estimated from the ellipticity pattern. The median FWHM and the ellipticity 
over the field is 0”.58 and 4.6 %, respectively. 

The median stellar image (PSF) size of each exposure yar- 
ied between 0”.52 and 0”.63 and the oyerall ayerage was 
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0”.58. Fig.[TJa) shows the map of the PSF size over the field 
of view (FOV). The non-uniformity is visible and we observe 
a hump of 0”.65 in lower left bottom and minimum of 0”.51 
at slightly upper-right from the center. This is not expected 
from the design and the tolerance analysis and suggests that 
at this early stage the collimation of the optical system was 
insufficiently accurate. 

To investigate this, we evaluate spatial variation of the el- 
lipticities of stars, {ei,e 2 }, which is defined as 

^ = {el,e 2 } = {/ii-/22,2 /i2}/ [in+hi + '^.J 

\ /(I) 

where Ijj is Gaussian-weighted quadrupole moment of the 
surface brightness distribution. 

Fig. Eb) shows a whisker plot where we present the posi¬ 
tion angle of the ellipse with a bar whose length is propor¬ 
tional to the ellipticity. We find a typical pattern of astigma¬ 
tism generated by the tilt of the camera optics with respect to 
the optical axis of the primary mirror. In this case, the cam¬ 
era is tilted around the axis shown in the figure as a dashed 
line. From the start to the end of this observation, the instru¬ 
ment rotator (InR) rotated by 68 degrees. The fix pattern seen 
on the FWHM and whisker plot co-rotated as the InR rotated 
which supports of our interpretation of the non-uniformity of 
the PSF. 

Recall that these test images were taken in the very first 
engineering run of the full FOV, we were trying to establish 
the way to measure the collimation error. Thus these particu¬ 
lar images were taken under rather poor alignment conditions. 
The median ellipticity over the FOV is roughly 5 % (Fig.Eb)) 
which is slightly larger tha n typical 3^4 % raw ell ipticity seen 
on Suprime-Cam images. (iMivazaki et al.l (^2007^ Fig.l). It is 
hard to determine the inclination angle from these data set 
alone but it is roughly a few arcmin level. By the time of 
this writing, we have established a way to measure the tilt and 
have implementered an auto-focus system, and we no lo nger 
see non-uniform patterns like this (IMivazaki et al.ll2E5h . 

2.2. Image Reduction and the galaxy catalog 

HSC has 116 2kx4k CCDs in total: 104 for imaging, in 
the central part of the field, and at the edges 4 CCDs for auto 
guiding and 8 CCDs for focus monitoring. We use the auto- 
guiding for the observations reported here. Each CCD has 
2048x4176 physical pixels (15/im: 0”.169 at the field cen¬ 
ter). The imaging area of each device consists of 4 512x4176 
stripes; each stripe is read through an independent output am¬ 
plifier located on one side of one of the shorter edges. Each 
segment has an 8 pixel wide pre-scan and we add 16 pixels 
of over-scan along the serial register and 16 lines of over-scan 
along the parallel register. 

We use the serial over-scan region for the zero reference of 
each amplifier (bias level). The median value of the pixel data 
on one row is used for the bias level for pixels on that row 
and subtracted. Non-physical reference pixels are trimmed 
from each segment and the segments are tiled together to re¬ 
construct the original 2048x4176 image. We then divide the 
image by a flat field image that is generated from an aver¬ 
age of dome Hatfield images (typically 10 images). Sky level 
is evaluated on by medianing on a mesh whose element size 
is 256x256 pixels and fitted with 2-D polynomial. This is 
subtracted from the image. This reduced CCD image file 
is one unit of the data source for the shape measurement 
(Sectior l231 l. 

The basic image analysis pipeline is developed through 


a collaboration of multiple institutions including Princeton 
University, Kavli IPMU and NAOJ. It is based on a cus¬ 
tomized version of the “LSST-stack” (llvezic et al. 1 120081 : 
lAxelrod et al.ll201(ih . a software suite being developed for the 
LSST project. Object detection is made on each CCD im¬ 
age. The PSP on a CCD is modeled as a function of the CCD 
pixel coordinates, and can be reconstructed accurately at any 
pixel position. This is used in the subsequent morphological 
classification process which tags star/galaxy flags as well as 
flagging cosmic rays. The flags are stored on the mask layer 
in an output FITS file. 

The star catalog is correlated with an external catalog 
(SDSS-DR8) for each CCD to perform photometric and as¬ 
trometric calibr ation, by using a c ross-matching alogorithm of 
astrometry.net (iLang et alJl2010l) . The SDSS magnitudes are 
transformed into the very similar HSC bands for photomteric 
zeropoints, and the astrometry.net is applied to the SDSS co¬ 
ordinates to determine the world coordinated system of each 
CCD. 
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Fig. 2.— The angular displacement of the control stars among exposures 
in the mosaic-stack process. The position averaged over the exposures is the 
reference (a) whereas the position in the external catalog is the reference (b). 

Because the number of stars on each CCD matched with 
astrometric catalog is limited (< 50), the accuracy of the as¬ 
trometric solution is not sufficient for the later image stacking 
process. In order to minimize the mosaic-stacking error, the 
residual of i-th control stars on a e-th exposure from a ref¬ 
erence frame, (Ajc/g, Ay,e), is parameterized as a polynomial 
function of the field position, (x,y), as 

N l N I 

= E E ^3'- = E E (2) 

/=0 m =0 /=0 m =0 

This i s a ’Jelly CCD’ approach adopted in IMivazaki et alJ 
(|2003), which is is particularly important for HSC where we 
have significant non-axisymmetric distortion pattern caused 
by the larger index mismatch o f the glass used for the at¬ 
mospheric dispersion corrector (IMivazaki et alj|2015l) . The 
implementation on the HSC pipeline is more sophisticated. 
Rather than dealing with the polynomial as a function of the 
pixel coordinate as above, the SIP convention (TAN-SIP) of 
World Coordinate System (WCS) is employed to allow the 
Jelly warp of the CCDs. A point on Fig.|2](a) shows the dis¬ 
placement of the astrometric position of the the control star 
with respect to the position on the first exposure, thus the scat¬ 
ter presents a measure of the mosaic-stacking error. This in¬ 
ternal error is as low as 6 to 7 mili-arcsec (mas). A point on 
Fig. |2] (b) shows the displacement from the external catalog 
and the scatter is larger (abut 60 mas) which we suppose is 
partly due to the astrometric error of the external catalog it¬ 
self. The peak position of the distribution in (b) has no offset 
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from the origin which indicates that there is no systematic er¬ 
ror on astrometric solution. The composition of the internal 
(a) and external (b) error gives a estimate of the total astro- 
metric error of ^ 20 mas in this analysis. Once the global 
solution over the multi exposures are obtained, the WCS of 
each CCD image is refreshed for further usage in the lensing 
analysis. 

Handling a entire mosaic-stacked image all at once is usu¬ 
ally not practical. We divide the image into a set of patches 
whose size is 4200 x 4200 pixel. In order to take care of ob¬ 
jects fallen on the boundary of the patches, we implement 100 
pixel overlap between neighboring patches. The patch is a 
unit of the following mosaic-stack operation and the object 
detection on the co-added image. We then combine the ob¬ 
ject catalogs of all paces into a single catalog by removing 
the objects detected multiple times within the overlapped re¬ 
gion. From the combined catalog, we extract the coordinate 
of moderately bright stars (22.0 < HSC-i < 24.5) for star cat¬ 
alog and the coordinate and the magnitude of galaxies. The 
star and the galaxy catalog as well as the image files with re¬ 
fined WCS of each CCDs (before the stack) are handed over 
to the next stage where we carry out the weak lensing shear 
estimate. 

2.3. Galaxy shape measurements and the shear estimates 

In the Sup rime-Cam weak lensing survey (iMivazaki et alJ 
[ 2002 i,[ 2003 ), we employed rather traditional method where 
we measured the galaxy shapes on the fully reduced mosaic- 
stacked CCD image. The PSF is not precisely round, nor is 
it of uniform size over the whole field. It is evaluated from 
star images over the whole field. (EqnlT]i. Because the PSF 
varies over the field of view (Fig. [Hb)), the PSF ellipticities 
are represented as a polynomial function of the field position. 
Analysis on the mosaic-stacked image requires that the PSF 
pattern does not change much during the series of exposures 
and that the CCD boundary does not cause a discontinuity of 
the PSF variation. Otherwise, the PSF becomes too compli¬ 
cated to represent as a simple polynomial function of the field 
position, and could cause systematic error in the galaxy shape 
measurement. In the case of the Suprime-Cam survey, we 
determined that the conditions were mostly met because one 
field was observed in a relatively short time span (40 minutes) 
with a small (^2 arcmin) dithering offset and we had a small 
CCD mosaic. However, in the planned HSC survey one field 
observation is split into several semesters to search for vari¬ 
ables and the circular FOV requires relatively large (nearly a 
half of the FOV) dithering step. It is therefore not expected 
that this simple requirement will be met. Instead, the galaxy 
shapes will be evaluated on the pre-stacked image and then 
combined after that to reduce the statistical error. When the 
survey is underway, we will develop code to determine shapes 
and fluxes simultaneously from the constituent images, which 
is arguably the most statistically efficient use of the data. 

For this early data, we have adopted Bayesian galaxy 
shape measurements implemented in the “/enifit” algorithm 
(iMiller et al.1120071: [Kitting et al.ll200i in this work. The 
PSF employed in the lensfit is not a simple elliptical but more 
empirical 32x32 pixel postage stamp image. Each pixel value 
of the postage stamp is fitted independently on each exposure 
to a 2D-polynomial function of the sky-coordinate, lensfit 
also allows varying low order coefficients between CCDs to 
further minimize the residual on each CCD. 

Galaxy shape measurement is made on the 40x40 pixel 
postage stamp image of each galaxy. Erom a galaxy coordi¬ 


nate on the given catalog, the postage stamp image is trimmed 
from the un-warped original CCD image by using the WCS 
information. Usually, the galaxy is imaged on multiple expo¬ 
sures and multiple postage stamps are created for one galaxy. 
lensftl has two galaxy model components; de Vaucouleurs 
profile bulge and exponential disk. The center of the two 
models are assumed to be aligned. The ellipticity, the 1-D 
size, the normalization and the bulge-to-disk ratio are the pa¬ 
rameters which model a galaxy. The postage stamp image of 
the model galaxy is convolved with the estimated PSE image 
at the galaxy position and the convolved image is compared 
with observed galaxy image to calculate the residual. By min¬ 
imizing the residual, the best fit galaxy model and the ellip¬ 
ticity is estimated. In the fitting process, multiple exposures 
data are simultaneously fitted. Note that the galaxy position 
is a free parameter as well, and is marginalized by integrat¬ 
ing under the cross-correlation function of the galaxy model 
and the data. Therefore, the astrometric error in the mosaic- 
stack process has limited impact on the error in the shape 
measurement. Note also that the observed image is not re¬ 
sampled nor warped onto the new coordinate in this process, 
w hich avoids the introd uction of correlated noise mentioned 
in iHamana et alJ (l2008h . 

Eor PSE modeling, we employ 2nd order polynominals for 
the exposure wide fit and 1 st order polynominal for the CCD 
term. Eig.[2(a) shows the ellipticity whisker plot of one expo¬ 
sure that is calculated from the PSE model postage stamp and 
averaged in the grid to visually compare with the observed 
PSF shown in Fig. [T] Fig. [2 (b) shows the difference of the 
observed and modeled ellipticity. The residual is sufficiently 
small as (Te < 1 %. 



Fig. 3.— Ellipticity whisker plot of modeled PSF (a) and the residual be¬ 
tween the model and the observed data. The sigma of the scatter is 0.7 % and 
0.6%, respectively in each component of the ellipticity. 


In the following sections, we adopt galaxies brigher than 
i = 24.5 for the mass map reconstruction. The rms sigma of 
the ellipticity of adopted galaxies is 0.41 and the galaxy num¬ 
ber density is 20.9 /arcmin^. In order to estimate the shear 
bias, we calculate mean ellipticities of galaxies used for the 
lensing analysis. The result is ((ei), (^ 2 )) = (0.005,-0.004), 
which is sufficiently small for this study. However, the bias is 
not completely negligible for cosmic shear studies and further 
investigations about the origin will be necessary in the future. 

2.4. Weak Lensing Convergence S/N map 

The dimensionless surface mass density, the convergence 
k(6) is evaluated from the tangential shear as 

K(0) = Jd^y,(cf>: 0 )Q(\cl)\) (3) 

where 7 ,(^ : 0) is the tangential component of the shear at 
postion 4> relative to the point 0 and Q is the filter function. 
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We adopt as a filter the truncated Gaussian 


Qc(d) = 


tt9^ 



(4) 


for 9 < 9o and Qc = 0 elsewhere (iHamana et alJlMT^ . We 
employ 6>c = 1 arcmin and 0^ = 15 arcmin, respectively. 

We adopt 9x9 arcsec grid cells and then calculate the con¬ 
vergence, k{6), on each grid using Eqn.© to obtain the k 
map. In order to estimate the noise of the k map, we random¬ 
ized the orientations of the galaxies in the catalog and created 
a Knoise map. We repeated this randomization 100 times and 
computed the rms value at each grid point. Assuming the k er¬ 
ror distribution is Gaussian, this rms represents the l-cr noise 
level, and thus the measured signal divided by the rms gives 
the signal-to-noise ratio (S/N) of the convergence map at that 
point. 

Fig. |4| is the convergence S/N map so obtained. 
IWittman et al.l(l2006b showed a convergence map based on the 
DLS R-band imaging. We find general agreement with that 
map although the resolution of the our map is higher . Note , 
however, that the convergence map of IWittman et al.l (l2006h 
was generated in the middle of the survey and it did not em¬ 
ploy the full depth image of DLS. They identifed two shear se¬ 
lected clusters on the map: DLSCL J0920.1H-3029 (140.033, 
30.498) and DLSCL J0916.0H-2931 (139.000, 29.526) and we 
see the corresponding peaks on our map. 

The former is a complex Abell 781 cluster region where 
at least four clusters have been identified so fa r from X-ray 
data and spectroscopic follow-up observations (ISehgal et al.l 
[200l . In Fig. |5] we showed a close up view of the 
region. Crosses indicate the center of X-ray emissions 
(XMM) of the four clusters which is called “West”(z=0.4273), 
“Main”(z=0.302), “Middle”(z=0.291) and “East”(z=0.4265) 
respectively from west to east. We detect clear corespon¬ 
dent peaks in the latter three clusters. This demonstrates that 
the angular resolution of the weak lensing convergence map 
matches well with that of the XMM X-ray image. 

We see two separate peaks on the convergence map at 
the “Main” cluster region and the mid-points of two peaks 
roughly coincides with the location of the X-ray emission cen¬ 
ter. The X-ray image is mostly round and no correspond- 
ing structure that i s found on the lensin g map is observed 
(ISehgal et al.ll200^ . IWittman et al.l (l2014t) presented the con¬ 
vergence map on this A781 region based on the full depth 
DLS imaging but no structure is seen in the “Main” cluster 
in that map. The structure seen in Fig. |5] survived several 
reality checks in our analysis (change of the magnitude cut 
of galaxies and bootstrap re-sampling of galaxies) but obvi¬ 
ously further observations are necessary to confirm the struc¬ 
ture. If it is confirmed, it will provide another laboratory for 
testing the natur e of dark mat t er foll owing the famous “bullet 
cluster”. In fact, ISehgal et al.l (120081) reported s mall extended 
sub-st ructure to the south-west of X-ray peak. iVenturi et al.l 
(1201 ll) discovered possible radio relic in their deep radio map. 
These are indicative of a merger and already suggested that 
the “Main” is not simple system. 

No significant peak is, however, found at the “West” po¬ 
sition in Fig. |5] In fact, the lensing signal of the X-ray 
emitting cluster Abell 781 “West” has been a matter of de¬ 
bate. Motivated by the less significan t detection on the orig- 
inal DLS map (IWittman et al.ll200^ . ICook & DelT Antoniol 
(1201 2|) examined the region based on images independently 
obtained by the Orthogonal Parallel Transfer Imaging Cam¬ 


era on the WIYN 3.5 m telescope and Suprime-Cam. They 
concluded that no significant signal was observed on either 
of the convergence maps and argued that this example could 
pose a challenge to the usefulness of weak lensing for the 
calibration of lensing mass against other obs ervable proper¬ 
ties of clusters in observational cosmology. IWittman et all 
(1201 4l) claimed, however, that the strong “joint constraint” 
of ICook & Dell’ Antoniol (1201 2l) was an artifact of incorrectly 
multiplying p-values. They re-visited the DLS data and tried 
to eliminate the contamination of foreground galaxies by 
adopting photometric redshift (photo-z) probability density 
weighting in the reconstruction of the convergence map. They 
suggested that the peak on the revised convergence map is still 
less significant (2.7 sigma level) but this does not strongly ex¬ 
clude the X-ray mass. We estimate the WL mass at X-ray 
position of the “West” cluster and the result i s shown in Ta¬ 
ble It is lower than what was obtained by IWittman et alJ 
(120141) but 1-sigma error overla ps each other. Theref ore, our 
result supports the summary by I Wittman et al.l (l2014h saying 
that the mass of “West” cluster (M 200 ) ranges 1-3 x 
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^^2000 

Fig. 5.— Close up view of the convergence map (Fig.|4j at the Abell 781 
multi-cluster region. Four cross marks are superimposed that indicate the 
location of the X-ray clusters observed by XMM: from west to east “West”, 
“Main”, “Middle” and “East”. 


2.5. Cluster Searches based on the Multi-Color Catalog 

The red symbols on Fig. |4] show the location of clusters 
of galaxies registered on NASA/IPAC Extragalactic database 
(NED). It is clear that the symbols tend to exist in the col¬ 
ored area (stronger lensing signal) and we see the general 
corespondence of the weak-lensing peaks and the cluster posi¬ 
tions. Among them, the star an d diamond symbols are SHELS 
clusters which are identified bv iGeller et al.l (120101) employing 
uniform spectroscopic observati on of a magn i tude-l imited (R 
< 20.6) sample using HectoSpec. iGeller etaH (120101) matched 
their clusters with the DLS lensing peaks; the star symbols 
show the matched clusters and the diamonds, un-matched. 

In order to make an independent comparison between the 
light and mass on this region, we search for clusters using 
the DLS public photometric data with a method, the “Cluster 
finding algorithm based on Multiband Identification of Red- 
sequence gAlaxies" (CAMIRA), developed by Oguri (2014). 
CAMIRA makes use of the stellar population synthesis (SPS) 
models of Bruzual & Chariot (2003) to compute SEDs of red- 
sequence galaxies, estimates the likelihood of being cluster 
member galaxies for each redshift using of the SED fitting, 
constructs a three-dimensional richness map using a compen¬ 
sated spatial filter, and identifies cluster candidates from peaks 
of the richness map. For each cluster candidate, the bright- 
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Fig. 4.— Weak lensing convergence S/N map reconstructed from the shear catalog. The smoothing radius 9g is 1 arcmin and galaxies brighter than 24.5 mag of 
HSC-i band are employed for the reconstruction. Red markers show the locations of clusters of galaxies registered on NASA/IPAC Extragalactic Dat abase where 
the obj ect type keyword ’GCluster’ was used to look up. The area below the declination of 31 degree is fully covered by Deep Lens Survey where [teller et al.1 
120101) carried out the spect roscopic campaign . Both Diamond and star markers show the location of SHELS clusters which match and un-match with he DLS 
lensing peaks, respectively IGeller et al.l2010() . 



est cluster galaxy candidate is identified based on stellar mass 
and location. Readers are referred to Oguri (2014) for more 
details of the algorithm and its performance studied in com¬ 
parison with X-ray and gravitational lensing data. 

We apply the CAMIRA algorithm to the DLS RV^z-band 
data. We apply a magnitude cut of R < 24.5 to exclude galax¬ 
ies with large photometric errors. In Oguri (2014), a number 
of spectroscopic galaxies in SDSS have been used to calibrate 
the SPS model. We adopt this SDSS calibration result, but we 


also add a constant 0.02 mag error quadratically to the model 
scatter, in order to accommodate the systematic offset of mag¬ 
nitude zeropoints and the difference in magnitude measure¬ 
ments between DLS and SDSS. We confirm that this SDSS- 
calibrated SPS model provides reasonable values when fit¬ 
ted to spectroscopic SDSS red-sequence galaxies in the DLS 
field. The cluster catalog contains, for each cluster, the clus¬ 
ter center based on the brightest cluster galaxy identification, 
the photometric redshift, and the richness. The redshft and 
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richness ranges are restricted to 0.1 < z < 0.8 and > 10, re¬ 
spectively. In addition, in this paper we also compute the total 
stellar mass by summing up weighted stellar mass estimates 
of individual galaxies; the weights are the “weight factors” 
Wmem^ which resemble the membership probability of each 
galaxy (see Oguri 2014). Stellar mass estimates of individ¬ 
ual galaxies are obtained from the SPS model fitting in which 
we assume a Salpeter initial mass function. 

We estimate the stellar mass only using optical data set 
which could cause the significant error in the estimation be¬ 
cause of the limited band-pass. When we argue the cluster 
members in this work, however, we mainly deal with passvie 
galaxies whose stellar mass is expected to be easier to esti- 
mate compared wit h that o f other general galaxies. In fact, 
lAnnunziatella et all ( 1201 4|) argued that the stellar mass of 
passive galaxies estimated from optical data agree well with 
what is estimated from optical and NIR data within 25 % if 
they adopt the template of passive galaxies in the calculation. 
Therefore we expect that our stellar mass estimates are not 
significantly biased by the absence of NIR data. 

In Fig. |2l we show the positions of resultant optically se¬ 
lected clusters as open circles where the center of the clusters 
is defined as the position of the bright cluster galaxy (BCG). 
The redshift is color-coded as is presented on the side bar. The 
diameter of the circles is 3 arcmin. 

3. RESULTS AND DISCUSSIONS 

3.1. Correlation of the Peaks on the Convergence Map and 
the Optically Selected Clusters 

We have searched for peaks in the convergence map; the 
results are indicated on Fig. [T] as filled triangles (significant 
peaks of S/N > 4.5) and filled squares (moderate peaks of 3.7 
< S/N < 4.5). We then match the peaks with the centers of 
the optically selected clusters and count the coincidences as 
we vary the match tolerance. As is shown in Fig. |6] when 
we increase the tolerance from zero, the number of matches 
rises rapidly and reaches a plateau at about 1.5 arcmin. The 
displacement of the BCG from the cluster center (as defined 
by the X-ray emission center) can be as large as 0.5 Mpc/h 
(IOgurill2014l) . which corresponds to ^ 1.5 arcmin at the red- 
shift of 0.7; the highest redshift in our sample. We would 
therefore expect to have to use a match tolerance of this order 
to recover real matches, and this is what we see; at approxi¬ 
mately this tolerance we have recovered all the real matches 
and the number should plateau, as it does. The number of 
matches increases again slowly beyond ^ 4 arcmin which we 
understand as accidental coincidences. 

Tahiti] shows the list of our peaks sorted by the S/N of the 
convergence map. Among nine significant peaks of S/N >4.5, 
five peaks (Peak ID 1, 3, 5, 6, 7) do not have a correspond¬ 
ing optically selected cluster using the CAMIRA algorithm 
described in the previous section. We carefully examine each 
case here. 

To complement our CAMIRA catalog, we looked for the 
clusters of galaxies in the NASA/NED database; the matched 
clusters are shown at the last column of the table where the 
tolerance is set at 1.5 arcmin. Peak ID 3 and 7 do have coun¬ 
terparts on NED (SHELS J0920.4h-3030;z=0.3004 and WHL 
J092104.1-i-303424;z=0.2758, respectively). These are proxy 
of the prominent A781 Main cluster (ID 0) and have similar 
redshift with the main cluster. In the CAMIRA algorithm, we 
eliminate the member galaxies of detected clusters to avoid 
double-counting. CAMIRA also uses a compensated spatial 


filter, which suppresses detection of clusters near very mas¬ 
sive clusters. These effects would explain why the CAMIRA 
catalog failed to detect clusters at the position of shear peak 
ID 3 and 7. The situation might be improved by modifying 
the form of the spatial filter for optical cluster finding but we 
leave this fine-tuning for the future work. 

Peak ID 5 matched with the original DLS shear s elected 
cluster; DLSCL J0916.0H-2931 (IWittman et al.ll200^ which 
is another complex system in this region. They reported that 
there are three associated X-ray peaks along north-south line 
and the north and the south peak were confirmed spectroscop¬ 
ically as clusters at a redshift of 0.53. We, in fact, note that 
below the richness threshold of 10 adopted here, there is a op¬ 
tically selected cluster at Photo-z of 0.542 (richness = 8.05) 
1.5 arcmin north of the Peak ID 5. Th e central X-ray pea k is 
later confirmed as a cluster at z=0.163 (iGeller et alJl2010h . 

Peak ID 6 is supposed be a possible sub-structure of the 
peak ID 0 because the X-ray emission peaks at just between 
Peak 0 and Peak 6 (Fig.|5]l. 

We therefore conclude that except peak ID 1 where no deep 
multi-color data is ayailable because it is outside of DLS field, 
all the other yery significant peaks of S/N > 4.5 are generated 
by physical entities. 

At a somewhat lower S/N leyel, iMiyazaki et al.l (120071) 
identified 17 peaks with S/N oyer 3.7 in a 2.8 deg^ region 
in XMM-LSS field, and found that nearly 80 percent of the 
peaks haye physical counterparts. On the other hand we haye 
26 peaks on the DLS oyerlapped 2 deg^ region and only 50 
percent of the moderately significant (S/N > 3.7) peaks haye 
identified physical counterparts. The discrepancy might be 
partially explained by the different noise leyel caused by the 
conseryatiye magnitude cut adopted in this work which in turn 
resulted in a smaller number density of weak lensing galax¬ 
ies, which in turn raises the noise leyel on the conyergence 
map. When we raise the threshold from 3.7 to 4.5 we see bet¬ 
ter match on the DLS field; nine peaks out of ten haye coun¬ 
terparts. Note that number density of the matched peaks are 
quite similar on XMM-LSS (4.3 peaks/deg^ for for S/N 3.7) 
and DLS (4.0 peaks/deg^ for S/N > 4.5). 



Fig. 6.— Number of shear selected clusters matched with optically se¬ 
lected clusters as a function of the tolerance of the angulai* distance used 
for the match. The the number of matches increases with the tolerance and 
reaches plateau around 1.5 arcmin and then gradually increases again, per¬ 
haps because of accidental matches. We set the tolerance at 2 arcmin (9 
matches). 
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TABLE 2 

Shear selected cluster samples generated by matching with the CAMIRA catalog and external catalogs. 


ID 

S/N 

RA2000 

DEC2000 

Photo-z 

log(M,) 

Richness 

Match 

0 

6.0 

140.082 

30.4898 

0.2807 

13.046 

96.872 

SHELS J0920.9+3029 (z=0.2915 A781 Main) 

1 

5.9 

138.998 

31.2972 




Out of DLS Field 

2 

5.7 

139.043 

30.4525 

0.6341 

12.594 

34.359 

Rank 2"* 

3 

5.6 

140.215 

30.4691 

0.3004' 



SHELS J0920.4+3030 (z=0.3004 A781 Middle) 

4 

5.4 

139.064 

29.8176 

0.514 

12.738 

38.957 

SHELS J0916.2+2949 (z=0.5343) 

5 

4.9 

138.993 

29.5623 

0.53F 



DLSCL J0916.0+2931 Rank P 

6 

4.9 

140.099 

30.5372 

0.2807 



Possible substructure of Peak ID 0 

7 

4.9 

140.256 

30.5689 

0.2758" 



WHL J092104.1 +303424 

8 

4.8 

138.982 

30.0472 

0.5204 

12.255 

12.372 

Rank 3^ 

10 

4.4 

140.299 

30.4687 

0.4022 

12.685 

41.929 

SHELS J0921.2+3028 (z=0.4265 A781 East) 

20 

3.9 

139.270 

30.0231 

0.3014 

12.583 

32.332 

WHL J091705.9+300118 (Photo-z=0.3285) Rank 0"* 

21 

3.9 

139.789 

29.5207 

0.326 

12.276 

10.285 

WHL J091906.0+293119 (Photo-z=0.3576) 

22 

3.8 

139.648 

29.4809 

0.544 

12.471 

11.963 


27 

3.7 

140.224 

29.6816 

0.274 

12.494 

26.238 

SHELS J0921.0+2942 (z=0.2964) 


^Geller et al.( 2010) 
HWittman et i 


‘ lUtsumi et al.l t20l^ 




A cluster (photo-z = 0.29) that is just below the richness 
threshold and has no associated peak is SHELS J0918.6+2953 
whose spectroscopically conhrmed redshift is 0.3 178. We no¬ 
tice th at this is one of the shear selected samples in iKubo et al.l 
(l2009l) (Rank=5 v = 3.9). Although we have positive con¬ 
vergence signal 2.5) here we see no strong peak (S/N > 
3.7). Also, we have a peak (S/N = 3.9) 5 arcmin north of 
SHELS J0918.6+2953. When we adopt a larger smoothing 
kernel {9g = 2 arcmin) on the convergence map, the positive 
signal is connected with the northern peak and results in more 
signihcant peak 3.7) which can be a counterpart of SHELS 
J0918.6+2953. In this paper, however, we will keep a single 
smoothing scale of 1 arcmin for easy comparison with theo¬ 
retical expectations. 

There is a cluster at photo-z of 0.5204 whose richness is 
low (12) but is mat ched with a simifica nt peak (ID = 8). This 
peak is reported in lUtsumi et al.l (I2014^ (Rank 3) as well us¬ 
ing a totally independent data set and data analysis pipeline; 
Suprime-Cam data analyzed by imcat. We therefore suppose 
that this is not a spurious peak caused by systematic errors. 
Although they observe a spatial concentration of eight galax¬ 
ies at redshift 0.537 (Az = 0.025), they did not identify the 
peak as a cluster becau se it did not meet the more stringent 
SHELS cluster criteria (iGeller et al.ll2010h . 


Fig. 7.— Peak location on the convergence map: filled triangle (S/N > 4.5) 
and filled square (3.7 < S/N < 4.5). Open circles sho w the locatio n of of 
optical clusters identified with the CAMIRA algorithm dOgurill2014ll whose 
richness is over 10 and open squares are clusters registered on NED (from 
west to east: CXOU J091554+293316 (z=0.184), SHELS J0920.9+3029 (z= 
0.291) and WHL J092104.1+303424 (z=0.2758). The diameter of the open 
circles is 3 arcmin. 

We now examine how the optically selected clusters match 
with the peaks on the convergence map in a different way. In 
Fig-El circles show the redshift and the richness of the clusters 
detected by CAMIRA algorithm. The clusters matched with 
the most significant peaks (S/N > 4.5) are marked with trian¬ 
gle symbols whereas the ones matched with moderate peaks 
(3.7 < S/N < 4.5) are marked with squares. It is encourag¬ 
ing to know that all luminous clusters samples (richness > 32) 
have corresponding convergence peaks. 



Photometric Redshift 


Fig. 8.— Photometric redshift versus richness plot of clusters detected by 
the CAMIRA algorithm (circles). Among them, clusters matched with the 
highest significance convergence peaks (weak leasing S/N > 4.5) are marked 
by triangles; lower significant ones (3.7 < S/N < 4.5) by squares. 
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3.2. Weak Leasing Mass Estimate 

The masses of the clusters are estimated from the tangen¬ 
tial alignment of the shear. The shear is azimuthally aver¬ 
aged over the successive annuli placed with a logarithmic in¬ 
terval in radius of 0.2. We fit the radial profile with a singular 
isothermal sphere (SIS) model, (surface density oc 0“'). In or¬ 
der to minimize dilution by the member galaxies of clusters 
(and any intrinsic alignment signal) we avoid the center; the 
fitting was made from the radius of 1.8 arcmin to 12 arcmin. 

The tangential shear, 77 -, induced by an SIS model is 


1 ^%s ^ 

Scr 2G r 


(5) 


where Y,cr, crsis and r is the critical surface density, the veloc¬ 
ity dispersion of the SIS and radial distance from the center; 


E 


cr 


Ds 

AttGDlDls^ 


r = DL0. 


(6) 


Here, D is angular diameter distance. By putting Eqn. ® into 
Eqn. (|5]l, we have 




(7) 


The mass of SIS within the radius, is given by 


M(< taJ = 



(8) 


The mass is also described by adopting the critical over¬ 
density parameter, A^,, with respect to the mean density, p(z), 
as 

Ttt t 

M(< taJ = —r\^p(z)Ac(z)/^M- (9) 

Eliminating ta^ from Eqn. ® and Eqn. (|9|l yields 


M(< =4.16 X l0^‘^MQh-Uasis/l000[km/s]f 

/ 1 500 ^ ( 10 ) 


where p(z) = (\+zf’LlMP% is adopted. We fit the data to obtain 
7 j-(0 = \ {deg\} and the mass is estimated by using Eqn. (|7]l 
and Eqn. (fTOl) where we emplo y the redshift distribu tion of 
the source galaxies, n(zs), from iLe Eevre et ^ (1201 3h rather 
than adopting the photo-z estimated only from optical BVRz 
data to avoid possible systematic effects from photo-z out¬ 
liers. Note that the weak lensing mass inferred from SIS 
model does not differ from what is obtained from NEW model 
in our fitting range. The difference depends on the mass but it 
is less than 10 % in most cases, which is rather smaller than 
the statistical error. 

Based on the compar ison of the CEHT M egaCam data with 
the GREAT simulation. iMiller et al.l (1201 3h argued that lensfit 
slightly underestimates the multiplicative factor, m, of the 
shear especially when the signal to noise ratio of the objects, 
V, is low; V <2Q whereas the underestimate becomes less than 
5 % when v > 30. Although we have not completed a com¬ 
prehensive comparison of HSC data with the simulation the 
behavior of lensfit on HSC data should not be significantly 
different from that of CEHT MegaCam because the image 
quality and pixel scale are similar. Therefore, we should have 
a level of roughly 5 % underestimate of shear at maximum 


TABLE 3 

Comparison of the mass estimate {Miqo/w'^Mq). The error is 
1(T. Note that West cluster mass is estimated at X-ray peak 
whereas all the other mass are estimated at the weak 

LENSING PEAK POSITION. 


A781 

This Work 

Wittman et al. (2014) 

Main 

Middle 

East 

7 

^■^- 1.6 

-1.2 

4 

‘t.O_l ^ 

6 7+1-4 

4.3fi:| 

2.8!!:^ 

West 

^-^-0.5 

_ 


which results in a roughly 7 % underestimate in the mass eval¬ 
uation. Since this is small compared with the statistical error, 
we will not deal with the systematic error explicitly in this 
paper. 

Table [3 shows the result of the comparison of our mass es- 
timate of A781 comp onents with the values in the literature 
(IWittman et aLllMTil) . Despite the totally independent obser¬ 
vation and the analysis, the agreement is encouraging and im¬ 
plies some progress in the convergence of weak lensing data 
analysis techniques. 



Fig. 9.— Richness versus Msoo^ of the shear selected samples that have 
optical counterparts detected by the CAMIRA algorithm (a). Redshift distri¬ 
bution of all identified shear selected clusters (b). 


Eig.|9](a) shows the relation between the richness and clus¬ 
ter virial mass, M 500 of the optically selected cluster samples 
that have counterparts with the convergence peaks. The cor¬ 
relation is clearly seen and the slope i s roughly 1 ^ 1.5, which 
agrees nicely with the slope found inlOgurj ([20l3) where the 
richness-mass relation was derived from stacked weak lensing 
analysis with the CEHTLenS shear catalogue. This further 
supports the reality of the cross-match of our lensing peaks 
and the optically selected clusters. 

InEig.|9](b), we show the redshift distribution of the iden¬ 
tified shear selected clusters on Tabl^ We see two spikes 
around z of 0.3 and 0.5 which clearly indicates that this nar¬ 
row 2.3 deg^ field is populated by large scale struc ture at these 
redshi fts; this has been mentioned already by iKubo et ^ 
(I 2 OO 9 II . We need a significantly wider field of view to over¬ 
come the local variance and to make cosmological arguments 
from the redshift distribution of clusters. 


3.3. Peak Count and Comparison with the Theoretical 
Estimates 
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The observed area which overlaps with DLS amounts to 
2.3 deg^, in which we found eight significant peaks whose 
S/N exceeds 4.5. Even if we drop Peak ID 6 from the list, 
which may be substruc ture within the A781 main cluster, 
seven peaks still remain. iHu & Kravtsov! d2003h estimated the 
cosmological variance in cluster samples and suggested that 
the variance exceeds the shot noise when the mass of clusters 
becomes less than ~ 3 x 10 '"^Mq with a weak dependence of 
the threshold mass on the survey volume. This is exactly the 
mass range that we are working on. Therefore, statistical ar¬ 
guments require comparison with cosmological simulations 
as we will see below. 
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Fig. 10. — Probability distribution of the number of the peaks under a given 
S/N threshold on a 2.3 deg^ wide field. Squares for S/N > 4.5 and triangles 
for S/N > 4.0. 

[Hamana et al.l (1201 2|) calculated the number of peaks on the 
weak lensing convergence map using a large set of gravita- 
tional lensing ra y-tracing simulations which are detailed in 
ISato et ail (l2008l) . Following the work, we made 1000 real¬ 
izations to evaluate the sample variance. In making the mock 
weak lensing convergence map, we added a random galaxy 
shape noise to the lensing shear data. The root-mean-square 
(RMS) value of the random galaxy shape noise was set so that 
the observed galaxy number density and RMS of galaxy ellip- 
ticities are recovered. We ado pted a fixed source redshift of 
Zj = 1.0. iLe Fevre et alJ (1201 3h estimated the redshift distribu¬ 
tion of magnitude limited samples taken from VIMOS VFT 
Deep Survey, and reported the mean redshift of (z) = 0.92 for 
a sample of 17.5 < Iab < 24 and (z) = 1.15 for a sample of 
17.5 < Iab < 24.75. Therefore, it would be appropriate to set 
Zj = 1 for our galaxy sample of Iab < 24.5. 

The expected peak count of S/N >4.5 is 0.61 on2.3deg^. In 
Fig . [TOl we show the probability distribution of the number of 
peaks under a given S/N threshold on a 2.3 deg^ wide field. As 
is shown in square symbols, the maximum number of peaks 
reached in the realization is 4 (9/1000) when the S/N = 4.5. 
This means that we have practically no chance to have seven 
or eight significant peaks on 2.3 deg^ field. Is this a challenge 
against the current CDM based cosmology? 

We note that the sensitivity of the number of peak to the S/N 
value is quite high, reflecting the steepness of mass function 
at the high mass end. So we experimentally lower the S/N 
down to 4 and examine the statistics. The mean number of 
peaks is 1.6 and the maximum number of peaks is eight in 
two realizations out of 1000 (0.2%, see triangles in Fig.fTOll. 
This still does not reconcile the gap between the observation 
and the prediction. 

In the mean time, iHamana et al.l (1201 2^ had adopted cos- 
mological simulatio n which used 3rd year WMAP result 
(ISnergel et alJl200^ . WMAP3 is known to have yielded a 


relat ively low q-r of 0 .76. If we adopt the recent Planck re¬ 
sult (lAde et al.ll2015h . erg = 0.83, the expected cluster count 
becomes a factor of 5.26 higher than WMAP3. This relaxes 
the tension dramatically and now what we observed is not ex¬ 
tremely unlikely; the chance to obtain more than 8 peaks of 
S/N > 4.5 is 3.7 %. One thing that we could suggest here is 
that our peak count strongly favours the recent Planck result. 

3.4. Stellar Mass Fraction in Clusters of Galaxies 



Fig. 11.— Fraction of stellar mass over the halo mass of the shear selected 
clusters (Circles). The dar k m atter virial halo mass, M^on, is estim a ted by 
putting Ac = 500 in Eqn. (To). Triangles are from (Gonzalez et al.l <20131) 
Table 6 and 7, and include the stellar mass associated with intracluster light 
(ICL). Solid line is the best fit to the triangle data. When they exclude the 
contribution of ICL outside of 50 kpc, the result is given by the dotted line. 
Typical error of 40 % is presented on the leftmost data point. Short-dash-line 
is the best tit power function for all the clusters and long-dash-line for the 
clusters excluding the most massive A781 Main. 

We now examine the ratio of stellar mass to halo mass, fs = 
Ms/Mhaio, of our samples. In the mass range that we probe, fs 
is reported to decrease as the halo mass increases, suggesting 
that star formation is less efficient in larger halos. This can 
be mostly explained by the inefflciency of th e cooling pro¬ 
cess inside larger halos. iGonzalez et alJ (I2013h presented one 
of the most recent results based on new halo mass estimates 
from XMM-Newton X-ray data. They claimed that the stellar 
baryon mass well compensates the shortage in the baryon bud¬ 
get and that the sum of the stellar baryon and the baryons in 
the form of gas almost reaches to the universal value estimated 
by WMAP and Planck. They also made a comparison of dif¬ 
feren t observational results (iLin et al.l 120031: IGonzalez et alJ 
l2007HAndreonll^l0tlFin et al.ll2012HFeauthaud et al.ll2012h . 

They found that the stellar m ass fractions reported i n other 
works are generally lower than IGonzalez et all (120131) . How¬ 
ever, they claimed that if the mass in the intra-cluster light 
(IGF) is considered, the discrepancy among the observations 
is minimized. 

Fig. [TT] shows the fs versus the halo mass of our samples 
(circles). The mass is estimated by weak lensing as explained 
previously and the mass range of the samples is wider than 
that of previous small samples. The stellar mass Ms shown 
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in Table 2 is the total stellar mass integrated by convolving 
spatial filter as shown in Fig. 2 of Oguri (2014). In order to 
estimate the stellar mass fraction accurately, we convert this 
stellar mass to the stellar mass within rsoo^ as follows. We 
assume that the stellar mass density profile follows an NFW 
profile with the concentration parameter as a function of halo 
mass of Duffy et al. (2008). For each halo with the mass 
4^500^, we derive a conversion factor from in Table 2 to the 
stellar mass within rsoo^ by convolving the projected NFW 
profile with the spatial filter of Oguri (2014) and estimating 
the ratio of the total mass with the spatial filter to Msoo^. We 
find that the conversion factor is ^ 1 for massive halos with 
4^500^ ^ 5 X IO'^Mq, and ~ 0.5 for less massive halos with 
4^500^ ~ 5 X 10^^4/q. We note that this correction also takes 
account of the 3D de-projection, i.e., properly removes mem¬ 
ber galaxies outside rsoo, projected along the line-of-sight. 
Thus our result shou ld be compared with t he so-called ‘3D’ 
stellar mass result in iGonzalez et al.l (1201 3h that incorporates 
the de-projection. The error in Mhaio estimate is roughly 30 %. 
We also expect a large scatter of ~ 30% in Mj for a given halo 
mass largely originating from the Poisson noise in the num¬ 
bers of cluster member and background galaxies. Therefore, 
we expect as large as 40 % level error in estimate of the in¬ 
dividual fs- The typical error bar is presented on the leftmost 
data point in Fig. |l2 _ 

Compared with lGonzalez et aT](l2013h results (triangles and 
the solid line in Fig. fTTIi . fs of this work is lower on aver¬ 
age over our mass range. We did not attemp t to include the 
contrib ution from the ICL but according to iGonzalez et al.l 
(1201 3h . the ICL contribution is ^ 25 % (dotted line in the fig¬ 
ure) which cannot totally explain the discrepancy. Our result 
agrees better with lower fs results es timated by Halo Occu¬ 
pation Distribution (HOD) method in lLeauthaud et akl (1201 2h 
and favors the argument that missing baryon problem has not 
yet been resolved in this mass range. We further note that 
IGonzalez et al.1 (120131) assumed a stellar mass to luminosity 
ratio that is on average slightly lower than the Salpeter ini¬ 
tial mass function adopted in this paper. Our stellar mass 
estimates are therefore larger than their estimates by ~ 10% 
which makes the discrepancy of the results larger by this fac¬ 
tor. 

The decrease rate of fs with increasing halo mass is an in¬ 
teresting observable because it reflects how the clusters are 
formed. Theoretical predictions and the recent N-body simu¬ 
lations coupled with semi-analytic models disfavor the steep 
slope in the framework of ACDM cosmological model. This 
is because the hierarchical clustering predicts that low mass 
halos with a certain fs would be assemble d together and be¬ 
come a larger halo with a comparable fs. iBalogh & McGeel 
(1201 Ol) suggested that the log-log slope of ~ -0.3 would be 
the upper limit to be consistent with their simulations in the 
ACDM cosmolo gy. The slopes that ha ve been observed and 
were compiled in iGonzalez et al.l (120131) span over -0.3~-0.6 
and the slope of -0.45 is presented by the author’s data which 
is slightly steeper than the theoretical prefere nce. In our case, 
the slo pe is slightly shallower than the data of iGonzalez et al.l 
(120131) : -0.32 when we exclude A781 Main (long-dash-line) 
and very shallow -0.05 (short-dash-line) when we consider all 
the clusters. It is, however, difficult to make really qualitative 
arguments here due to the large error for the limited number 
of clusters in this work. More data are certainly cried for and 
the on-going HSC legacy survey will provide ideal sources for 
the future studies. 


4. CONCLUSION 

We show the results of a weak-lensing cluster search on 2.3 
square degrees of HSC commissioning data in the Abell 781 
field with 1.6 hours of exposure. The data are excellent, with 
very good image quality; some slight astigmatism was found 
which was traced to a small miscollimation (since corrected) 
but was not difficult to correct because of the low order spatial 
variation. 

Clusters were searched for on a high resolution conver¬ 
gence map generated froma these data. We see very good 
agreement with the previous results in mass measurements in 
this field made by Wittman et al (2014) except for one peak, 
corresponding to Abell 781 West. Clusters of galaxie s were 
search ed for in optical data using CAMIRA algorithm (iQguril 
l2014h . This cluster list was compared with the locations of the 
peaks in the convergence map. There is only one significant 
peak which we cannot judge the reality because it is outside of 
DLS field where no multi-color data is available. All the other 
peaks of S/N > 4.5 are physically real. This demonstrates the 
reliability of the convergence map generated from these even 
early HSC data. The number of observed clusters at this level 
is significantly larger than the predicted average number in 
this field size expected from theoretical calculations based on 
WMAP-3 results. This result, however, is extremely sensi¬ 
tive to the value of ag in the theoretical predictions, and is 
not unlikely if we adopt the recent Planck cosmology with its 
somewhat higher value of trg = 0.83. 

The stellar mass fraction in our sample is systematically 
lower than one of the other most recent results and is more 
consistent with the earlier values estimated by use of the HOD 
statistical formalism. Our result thus favors the argument that 
baryons are still missing in this mass range. The decrease of 
the stellar mass fraction with increasing halo mass is slightly 
shallower than the previous work, and is more consistent with 
current (though still very uncertain) simulations. 

Because of both the limited sample size and the statistical 
errors in our masses, the results for the stellar mass fraction 
are not strictly conclusive, but are strongly suggestive. 

So even though the results for number density and stellar 
mass fraction from this very small sample are not conclusive, 
samples not enormously larger than this with this instrumen¬ 
tation will be sufficient to provide excellent results on these 
important quantities, and a large survey is underway which 
will provide these data. 

In this work, we were at least able to demonstrate that 
cluster identification, redshift estimates, and mass estimates 
can be obtained by multi-band optical imaging data with the 
newly developed Hyper Suprime-Cam camera through weak 
lensing and cluster finding techniques. HSC has uniquely 
combined features of wide field, large aperture and superb 
image quality and the data from the currently on-going HSC 
legacy survey is very promising for these and many other cos¬ 
mological investigations. 
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